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1 Introduction 

One of the peculiar features of the two dimensional four state Potts model 
is the presence of a marginal field which leads to universal multiplicative 
logarithmic corrections to the scaling laws. 

These corrections can be exactly evaluated [0, §] but, as often happens 
when dealing with marginal fields, they are accompanied by still large sub- 
leading non-universal contributions which completely mask the scaling be- 
haviour of the system (at least for those values of the correlation length which 
can be reached in standard simulations). 

This is different from the behaviour, for instance, of the 3d Ising model 
where even at moderate values of the correlation length, the non-universal 
corrections [§] give very small contributions; thus they can be safely taken 
into account by adding to the scaling functions only the first non-universal 
term. 

In the present case, instead, these corrections are so large that the re- 
liability of the fits, in which only the first subleading term is taken into 
account, becomes questionable. On the other hand, with current numerical 
precision, it is almost impossible to add further corrections without losing 
any predictive power in the fits. 

This makes the numerical study of the four state Potts model one of 
the most difficult tasks in the context of Monte Carlo simulations of two 
dimensional spin models. 

This problem was recently addressed by J.Salas and A.D.Sokal in @] by 
extending the RG analysis of [^ ^ up to third order in the fields. They 
succeeded in obtaining the universal leading corrections to the scaling, which 
turn out to be additive terms of the generic form log log /log. The improved 
scaling functions were then tested by looking at the critical finite size prop- 
erties of the model, and an improvement of the scaling behaviour of the 
data was observed which however turned out to be still affected by large 
non-universal 1/log terms. 

In this paper we return to this problem by looking at various universal 
amplitude combinations of the model |^. We find an improvement in the 
scaling behaviour of our data if the universal contributions evaluated in [Q 
are taken into account. 

However, as in |^, this is not enough to describe the data; non-universal 
corrections must be considered and final results crucially depend on the type 
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of terms included in the scaling functions. We consider this one of the most 
delicate aspects of this paper. For this reason we described, as precisely as 
possible, the procedure used to construct the scaling functions (see sect. 7) 
and included in the paper, besides final estimates for the amplitude ratios, 
also the direct results of the Monte Carlo simulation (see tables 1,2 and 3) so 
that the reader can use the data to study alternative scahng functions and 
possibly find a clever way to control the systematic errors involved in the 
truncation that we suggest. 

Fortunately for the present problem we have an independent way to test 
our results. In fact, thanks to a recent work by J.L.Cardy and G.Delfino P, 
precise estimates for amplitude combinations are now available in the con- 
tinuum limit. The relatively good agreement found between our estimates 
and those of make us confident of the reliability of our results and at the 
same time strongly supports the correctness of the form-factors derivation of 
ref. i. 

This paper is organized as follows. In sect. 2 we give a general intro- 
duction to the g-state Potts model and we summarize a few known facts 
concerning the phase diagram and its scaling limit description in the frame- 
work of conformal and perturbed conformal field theory. In sects. 3 and 4 we 
introduce and discuss the observables and the amplitude ratios in which we 
are interested. In sect. 5 and 6 we describe the Monte Carlo simulation and 
test the results by comparing them with existing low temperature series and 
by imposing the duality relations on the internal energy and the specific heat. 
Sect. 7 is devoted to the study of the scaling behaviour of the observables 
and to extract the best estimates for the amplitude combinations. Sect. 8 is 
devoted to a comparison of the results with those of ref. while in sect. 9 
we have collected some concluding remarks. 

2 The Model 

We study the four state Potts model in two dimensions on a simple square 
lattice. The action is given by 

Spotts = —(3 2^ 5s{x),s{y) , (1) 

<x,y> 

where the field variable s{x) takes the values 0, 1, 2, 3; x = (xq, Xi) labels the 
sites of the lattice and the notation < x,y > indicates that the sum is taken 



on nearest neighbor sites only. The 6 function is defined as usuah 6a,b = 1 if 
a = b and otherwise. The couphng (3 is related to the temperature in the 
standard way (3 = -gp. In the following we shall always consider lattices of 
equal extension L and periodic boundary conditions in both directions. 

Several results are known exactly for this model. The action eq. (|l]), is 
invariant under the permutation group 5*4. However in the low temperature 
phase this symmetry is spontaneously broken to S3. The two phases are 
related by duality and separated by a second order phase transition located 
at Pc = ^ = log(3) = 1.098612. The dual coupling (3 is related to the 
original coupling /3 by 

and the fixed point of this relation is the critical coupling /5c. 
It is useful to introduce the variables 

(Ta{x) = 6s(x),a - J , = 0,1,2,3 . (3) 

It is easy to see that < ctq, >= Va in the high temperature phase and 
that they all become different from zero in the low temperature phase. In 
particular for one of the four values of the spin (which we shall call in the 
following "majority spin") < a^ > > 0, while for the three other values we 
have < aa > < 0. 

It has been shown (cf. Q) that the partition function Z{T,q) for the 
g-state Potts models {q integer) on a square lattice A can be written as 

Z(T,g) = 5:(e^-l)V • (4) 

{G} 

The sum in eq. (^ runs over all the graphs G on A, C is the number of 
connected components (including isolated sites) in G, and A^ is the number 
of bounds on the lattice edges. For a better understanding of some of the 
peculiar features of the four-state Potts model it is convenient to consider 
the phase diagram of the whole family of models, defined for arbitrary q > 
0. Eq. (0) provides an expression suitable for extending the definition of 
Z(T, q) to non integer values of q. The g-state Potts model undergoes a 
phase transition at 

/3c = log(v/g + l) . (5) 



Below this temperature the system is in its Sg-broken symmetry phase whereas 
above it the system is fully disordered. The transition at T = Tc is first order 
for g > 4 but becomes second order for g < 4, in the latter case the model 
gets renormalized into a conformal field theory with central charge 



^=i-(7^ • («' 



where / is related to q by 



In the scaling limit, at rational points / = p/{p' — p), the thermal field e 



rescales with scaling dimension [T^ 



A.4(l + ^) . (8) 

hence it can be identified with the operator 02i in the Mp^ minimal confor- 
mal model. Notice also that eq. (^ shows that two square-root branch points 
at g = and q = 4 are present. At g = the thermal operator is marginal, 
in the analytically-continued second branch it becomes irrelevant and the 
critical point has moved into the antiferromagnetic region. The physics in 
this sector is certainly very interesting but it is slightly less relevant for our 
current interests. More related to the subject of this paper is, instead, the 
physical meaning of the second branch point at g = 4. Let us consider a 
further variant of the model in which vacancies are allowed, and correspond- 
ingly a chemical potential n is introduced. In the sector < g < 4 with 
H negative or sufficiently small, the additional dilution field turns out to be 
irrelevant and the system still undergoes a second order phase transition in 
the same universal class of the pure Potts model. Near the transition point 
the dilution field scales with a conformal dimension 

A/. = 2 + ^ , (9) 

and it can be identified with primary conformal operator </)3i. 

At g = 4, the dilution field (p^ becomes marginal; along the critical RG 
fiow its slow rate of disappearance cause now multiplicative logarithm cor- 
rections to the critical behaviour. From eq. (0) we also see that the net result 



of the entrance in the second branch consists of a negation of /. Hence at 
the same value of g, but on the second branch we now have a conformal field 
theory (CFT) with central charge 
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/(/ + i: 



with thermal and dilution fields with dimensions 



A, 



1 



/ + ! 



A, 



/ + ! 



(10) 



(11) 



These two fields can now be respectively identified with the relevant confor- 
mal operators 0i2 and 0i3. In conclusion, at fixed (? < 4, the phase diagram 
in the plane (yU, T) is as follows: if /i is negative or sufficiently small then 
at T = Tc{fi) the system undergoes a second order phase transition in the 
universal class of the pure Potts model, whereas when /x is large the tran- 
sition becomes of the first order. On the critical line (/i,Tc(yu)) the point 
marking the change of critical behavior is the tri critical point. From the 
scahng quantum field theory point of view, the picture looks also consistent 
with the one described above. First notice that perturbation of the confor- 



mal field theory with either the thermal or the dilution field is integrable []TI 



and that the associated quantum field theories have been the subject of very 
deep studies (see for example [0, |1^, |1^, |15|, |1^, |T^ ) . The thermal operator 



e drives the system into a massive ordered or disordered phase depending on 
the sign of the perturbing parameter |]15|, |1^, [l^ . The operator 0^ instead 
moves the tricritical model either into a massive phase (a line of first-order 
transitions) or into a critical massless phase. The IR fixed point of the latter 
is the Potts-model CFT and the two less irrelevant attracting operators are 
the fields 03i and TT O, O, HI . 



Quantum reductions of the Izergin-Korepin model at rational points, 
giving the S-matrix elements for the 0i2/02i perturbations, were first ob- 



tained by F.A.Smirnov [0. Subsequently using a somehow different kink 
basis, L.Chim and A.B.Zamolodchikov [|TB| defined alternative scattering el- 
ements for the model. The latter formulation, being more suitable for ana- 
lytic continuation at arbitrary values of q, has been used by J.L.Cardy and 
G.Delfino to make predictions about the values of some universal ampli- 
tudes. The method used is a variant of the form-factors approach to the 
correlation functions proposed in p!9 |. 



3 The observables 

3.1 Magnetization 

The magnetization of a given configuration is defined as: 

m = -J2(^c.n.{x) ' (12) 

where \^ = L^ is the volume of the lattice and am is the value of spin corre- 
sponding to the majority of the spins. However, in a finite volume at arbitrary 
finite low temperature the 5*4 symmetry of the model is not spontaneously 
broken. Practically this means that, in the simulation sample, configurations 
with all the four possible values of am appear with equal probability. In order 
to obtain a low temperature non-vanishing magnetization a magnetic field 
h that explicitly breaks the symmetry, must be coupled to the system. The 
thermodynamic limit at non zero h should be taken first, then the limit of 
vanishing magnetic field could be performed. However it is difficult to follow 
this route in a numerical study. An alternative, commonly adopted, approach 
is to identify am in each configuration by counting the spins belonging to the 
four possible values of a and then extracting the majority one. This pro- 
cedure works in a satisfactory way if the lattice size and coupling constants 
are such that the probability of finding interfaces among different vacua is 
negligible. We carefully chose our lattice sizes so as to satisfy this bound Q. 

In the following we shall assume for simplicity that am = is the value 
of the majority spin and shall denote the remaining three values with roman 
indices i,j, k, ... = 1, 2, 3. 

Close to criticality and at t = ^^^^ < 0, the magnetization scales as ||^ 



^Let us note, as a side remark, that this procedure, in the Ising case is equivalent to 
the choice 

< m > = hm < |m| > (m = — > Si) , 

L-->oo V ^ — ^ 

i 

where the s^'s are in this case Ising spins. 

The finite size behaviour of this observable was carefully studied in pQ] . It was shown 
that this choice converges to the infinite volume value better than any other existing 
proposal and that the asymptotic, infinite volume, value is reached for lattices of size 
L >~ 8f , where f denotes the correlation length. In our simulations we always used 
lattice sizes much larger than this threshold. 



<m>^ i? (-t)i2(-log(-t)) 8 



• 3 1og(-log(-t)) , 

16 -log(-t) ' 



log(-t). 
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3.2 Magnetic susceptibility 

The susceptibility 



X 



d < m > 
dH 



(14) 



gives the response of the magnetization to an external magnetic field and it 
can be expressed in terms of moments of the magnetization 



< m > 



< m>' 



x = v{< 

In the high temperature phase this means 

')^ = V<m^> =V<al> 



(15) 



(16) 



where a is any one of the four values (0, 1, 2, 3)^ 

In the broken symmetric phase, depending on the choice of coupling the 
external magnetic field to the majority spin or to one (or more) of the other 
values, two kinds of susceptibilities can be defined p2[ 



Xi 
Xt 



<o^o > 



< i(^^ - (^j? > (^ ^ j] 



longitudinal 

transverse 



(17) 
(18) 



In this paper we concentrate on the longitudinal susceptibility xi- 



^Notice however that, inspired by the analogy with the Ising model or by the embed- 
ding in the AT model, different definitions of the order parameter (hence of the mean 
magnetization in the symmetric phase) are possible. For instance 

< TOi >=< aa — ap > 

with a 7^ /3 or 

< 1712 >~< CTq — (T/3 + cr^ — (TS > 

with a ^ P j^ J ^ S. The corresponding susceptibilities are related to x by simple 
multiplicative constants. In comparing our results with those of pl| , H one must take into 
account this different normalization. 



Close to the critical temperature x and xi scale a^ 



r 9iog(-iog|t) . J ^ \ 

[ 8 -log|t| "^^logltlj 

r 9iog(-iog|t) ^q( ^ \] 


(t>0) 
(t<0) 


8 -log|t| Vogt) 



x~r+(t)-6(-iog|t|)4 

X;~r_(-t)-i(-log|t|)i 



3.3 Internal energy and specific heat 

The internal energy is defined as 



(19) 



(20) 



the specific heat 



C ^ ^^^ = 2V(<E'>-<E>' 
dp ^ 



(21) 



Duality relates both internal energy and specific heat in the low temperature 
phase to those in the high temperature phase. The relations are 

{l-e-^)E{P) = l-il-e-^)EiP) , (22) 

C{(3){l-e-^y+e-^{l-e-^)E{f3) = C0){l-e-^y+e-^{l~e-^)E0). (23) 
Close to the critical temperature we have 

3 log(- log |t| 



C~ A+ (t)-3(-log|t|)-i 1-^- 
C~A_ (-t)-i(-log|t|)-^ 



2 -logltl 



+ 01 



^ 3 1og(-log|t|) ^ ^ 1 



Tog t 



1 Y 


(t>0) 


\og\t\J_ 


1 \ 


(t<0) 


log|t|j_ 



(24) 



^Notice that there is a misprint in the analogue of this equation in ref. 



3.4 Second moment correlation length 

We consider the decay of so-called time-slice correlation functions. The mag- 
netization of a time-slice is given by 



Sa{Xo) = Y'^Cra{Xo,Xi) 



XX 



Let us define the correlation function 

Gafiij) = ^{{Sa{Xo)Sf3{Xo + t)) - (S'„(Xo)) {Sp{Xo))} 



(25) 



(26) 



Xo 



For any choice of the indices a and (3 in (^6]), the dominant large distance 
behaviour of G{t) is dominated by the lowest mass of the model: 



G{t) oc exp(-r/0 



(27) 



where ^ is the exponential correlation length and coincides with the inverse 
of the lowest mass of the model. However, at low temperature, the rich 
structure of the spectrum (cf. 0) can mask this asymptotic behaviour. For 
this reason, in the following, we shall concentrate on Goo which has the 
"neatest" asymptotic behaviour. This will also allow us to directly compare 
our results with those of [§. In the high temperature phase the lowest mass 
is well separated from all other excitations in the spectrum and extracting 
the exponential correlation length is much simpler. 

Close to criticality the behaviour of the correlation length is governed by 
the scaling laws 

^r^U{t)-H-\Og\t\f^ 

^^/_(_t)-|(_iog|t|)i 

^ — log \L\ \ log \i\ / 

(28) 
We are also interested in the second moment correlation length ^2, to evaluate 



^ 1 31og(-log|t|) 1 c.f Ml 
4 -logt| Vogt) 


(t>0) 


^ 1 31og(-log|t|) 1 ^/ 1 \1 
4 -logt| Vogt) 


(t<0) 



^2 we used the estimator [21 



e 



(2) = 



i-1 



2sin(7r/L) 



(29) 



where x is the susceptibihty and F is the Fourier transform of the correlation 
function at the smallest nonzero momentum (27r/L,0). 

Notice that the susceptibility can be rewritten as zero momentum Fourier 
transform of the correlation function, hence, in order to have a consistent 
definition, the same correlation function must be chosen in both x ^i-nd F. 
In the low temperature regime we are interested in setting the longitudinal 
susceptibility 

X, = EGo,o(x) (30) 



in eq. (|29D , hence, we shall study 

F = 5]e2--i/^Go,o(a:) • (31) 

X 

^2 is a very popular approximation for the exponential correlation length 
since, in Monte Carlo simulations, its numerical evaluation is much simpler 
than that of ^. 

Moreover it is the length scale which is directly observed in scattering ex- 
periments. It is important to stress that ^2 and ^ are not fully equivalent (cf. 
P ), even though their critical behaviours are the same up to a multiplicative 
factor. In particular, the ratio ^/^2 gives an idea of the density of the lowest 
states of the spectrum. If the lower excited states are well separated the 
ratio is almost one, whereas a significantly bigger ratio indicate a denser dis- 
tribution of states. Furthermore, different choices of the correlation function 
in eq.s (|30|j3T|) lead to different values of ^2, while the exponential correla- 



tion length is always the same. A careful study of these differences can give 
several information on the spectrum of the theory. We shall call the critical 
amplitudes /2,± for ^2 to distinguish them from /±. 

4 Amplitude ratios 

We are interested in the following amplitude ratios 
and the following amplitude combinations 

r r 
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which are scale invariant thanks to the (hyper) scahng relations among the 
critical exponents 

a . (34) 



a + 2/? + 7 = 


= 2 


du = 2- 


are also interested in the combinations 


R^ - ^-^^^ 


5 


R+ = ^A+h, 


^c - ^2 


R- - ^-^- 


? 


R^ = ^A.h, 


^^ - R2 



(35) 

(36) 

which have particularly interesting behaviours in the 4 state Potts model (see 
below). We shall neglect the amplitude ratio A^/A_ which is trivially 1 due 
to duality 0. 

5 The simulations 

We produced a standard cluster algorithm using both the Wolff single cluster 
update and the Swendson Wang cluster update. After prehminary tests, we 
used the latter algorithm for our high statistic simulations. 

To check our program we made comparisons of the MC results with the 
exact solution on a 3^ lattice and with the Salas and Sokal (ref.s |^l], ^) 
results at the critical point on a 16^ lattice, with a comparable statistics. 

We simulated the 4 state Potts model in the high and low temperature 
phases for 16 values of the couplings which were chosen exactly as dual pairs. 
This choice allowed us first, to perform a very stringent test on our estimates 
for the thermal observables which must be related by duality and second, 
to obtain a direct estimate of some amplitude ratios. The results for the 
observables in which we are interested are reported in tables 1,2 and 3. 

Lattice sizes were chosen large enough to make finite size effects negligi- 
ble within our statistical errors. After a preliminary test on the finite size 
behaviour of our observables we chose L > 10 ^ in the high temperature 
phase and L > 20 .^ in the low temperature phase, a default size L = 120 
was taken for small values of C,- In each simulation the number of measure- 
ments was 2 ■ 10^. Each measurement was separated from the next one by 



^In principle this result could be used to test our simulation. But such test is completely 
equivalent to the test of the duality relation eq. (E3) that we perform in tables 4 and 5. 
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two Swendsen-Wang updates. A standard jacknife procedure has been used 
to analyze statistical errors. 



p 


L 


^2nd 


E 


C 


X 


1.06722 


120 


6.44(5) 


0.64604(6) 


1.571(12) 


11.83(2) 


1.07722 


120 


8.46(4) 


0.66318(6) 


1.892(13) 


18.78(3) 


1.07972 


120 


9.25(3) 


0.66808(6) 


2.018(14) 


21.80(4) 


1.08222 


120 


10.26(3) 


0.67337(6) 


2.189(15) 


25.90(5) 


1.08472 


120 


11.60(3) 


0.67909(6) 


2.405(15) 


31.69(6) 


1.08722 


180 


13.37(5) 


0.68544(6) 


2.614(22) 


40.31(9) 


1.08972 


180 


15.94(5) 


0.69232(6) 


3.033(24) 


54.05(13) 


1.09222 


240 


20.33(7) 


0.70056(6) 


3.574(35) 


80.88(24) 



Table 1: Results in the high temperature phase. 



p 


L 


^2nd 


E 


C 


1.130500 


120 


2.85(5) 


0.850888(5) 


1.3718(10) 


1.120231 


120 


3.69(5) 


0.835061(7) 


1.7474(16) 


1.117684 


120 


4.00(4) 


0.830453(8) 


1.8851(17) 


1.115135 


120 


4.45(4) 


0.825436(8) 


2.0516(20) 


1.112597 


120 


5.01(4) 


0.820001(9) 


2.2604(24) 


1.110065 


120 


5.77(4) 


0.813951(10) 


2.5274(28) 


1.107540 


120 


6.98(4) 


0.807100(12) 


2.9048(40) 


1.105020 


180 


8.89(7) 


0.799095(9) 


3.4872(49) 



Table 2: Results in the low temperature phase: thermal observables and 
correlation lengths. 



6 Analysis of the results 
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p 


L 


m 


Xi 


1.130500 


120 


0.633288(7) 


1.2180(9) 


1.120231 


120 


0.612874(11) 


2.1072(23) 


1.117684 


120 


0.606331(13) 


2.5080(26) 


1.115135 


120 


0.598834(14) 


3.0572(36) 


1.112597 


120 


0.590219(16) 


3.8447(51) 


1.110065 


120 


0.579950(19) 


5.0315(72) 


1.107540 


120 


0.567272(26) 


7.053(14) 


1.105020 


180 


0.550686(23) 


10.987(19) 



Table 3: Results in the low temperature phase: magnetic observables. 



6.1 Energy and specific heat 



By using eq. (|2^) we have a non trivial test of our estimates both for the 
energy and for the specific heat. In table 4 we compare the results for the 
internal energy in the high temperature phase with those obtained using 
eq. ( p2D and the values measured with the dual coupling, at low temperature 
as input. A similar comparison for the specific heat can be found in table 5. 



6.2 Magnetization and Susceptibility 

In tables 6 and 7 we compare our results for the magnetization and the 
low temperature longitudinal susceptibility with the series of ref. p3|. Both 



for the magnetization and for the susceptibility we used the diagonal Fade' 
approximant. As expected, the agreement, which is rather good far from 
the critical point becomes worse and worse as (3c is approached. Notice 
however that we used the simplest possible resummation technique, more 
sophisticated approaches like the double biased IDA of ref. p4| could give 
better results and could also give a way to estimate the systematic errors 
involved in the truncation and resummation of the series (for an attempt 
in this direction in the case of the 3d Ising model see for instance 0)0. In 



^In comparing our values of the magnetization with those of ref. p3| one must notice 
that there is a factor 4/3 between the two definition of magnetization. On the contrary 
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/? 


HT 


LT+eq.(il) 


1.06722 

1.07722 
1.07972 
1.08222 
1.08472 
1.08722 
1.08972 
1.09222 


0.64604(6) 
0.66318(6) 
0.66808(6) 
0.67337(6) 
0.67909(6) 
0.68544(6) 
0.69232(6) 
0.70056(6) 


0.646061(5) 
0.663179(7) 
0.668074(8) 
0.673365(8) 
0.679055(9) 
0.685341(10) 
0.69240(1) 
0.70060(1) 



Table 4: Comparison between the internal energy measured in the simula- 
tion at high temperature (second column) and the values obtained, from the 
internal energy measured at low temperature, using the duality relation 
(third column). 



the case of the susceptibihty the discrepancy between the results from the 
series expansion and our Monte Carlo are larger and only the first value of 
beta agree within the errors. It is clear however that we are pushing the 
series to their limit of validity and in fact, looking at nondiagonal Pade' 
approximants, one sees very large fluctuations (much larger than in the case 
of magnetization) in the series estimates. 

7 Scaling behaviour 

Let us now address the problem of extracting the continuum limit values 
of the quantities discussed in the previous section. As mentioned in the 
introduction, due to the presence of large corrections to scaling terms, this 
requires a rather non-trivial analysis. We followed a three step procedure. 

1] As first test we tried to fit the data using only the dominant multi- 
plicative log correction keeping into account for C and x the possible 
existence of bulk constant terms. Hence a one parameter fit for ^2 smd 
m and a two parameter fit for C and x- Iii ^ill cases we found very 



there is complete agreement between the two definitions for the longitudinal susceptibility. 



14 



/? 


HT 


LT+eq.(|3D 


1.06722 


1.571(12) 


1.555(1) 


1.07722 


1.892(13) 


1.903(2) 


1.07972 


2.018(14) 


2.032(2) 


1.08222 


2.189(15) 


2.191(2) 


1.08472 


2.405(15) 


2.391(3) 


1.08722 


2.614(22) 


2.646(3) 


1.08972 


3.033(24) 


3.011(4) 


1.09222 


3.574(35) 


3.579(5) 



Table 5: Comparison between the specific heat measured in the simulation at 
high temperature (second column) and the values obtained, from the internal 
energy and the specific heat measured at low temperature, using the duality 
relations (22, '2^ (third column). 



high Xr and, even ehminating all the data except the two couplings 
nearest to f3c, it was impossible to reach a reasonable confidence level. 
This clearly indicates that additive corrections to the scaling cannot be 
neglected. 

2] The second step was to add the first non-universal correction, (that 
of the form 1/log ). At the same time we also included the universal 
corrections evaluated by Salas and Sokal, which are of the same order of 
magnitude and do not add degrees of freedom in the fit. The resulting 
fitting functions are (see sect. 3): 



x(t) = ao + r±|t| 



log|t|) 



1 + 



91og(-log|t|) ai 



m{t) 



f2,±\t\ 



5|t|T2| 



;-iog|ti)' 

log|t|)-5 



- log \t 

3 1og(-log|t|) 
4 -log|t| 

3 log(-log|t|) 



+ 



+ 



Oi 



log|t| 



16 — loffltl 



+ 



■log|t| 

ai 

-log|t| 



(37) 
(38) 

(39) 
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(3 


our MC 


series 


1.130500 
1.120231 
1.117684 
1.115135 
1.112597 
1.110065 
1.107540 
1.105020 


0.633288(7) 
0.612874(11) 
0.606331(13) 
0.598834(14) 
0.590219(16) 
0.579950(19) 
0.567272(26) 
0.550686(23) 


0.633275 
0.612863 
0.606310 
0.598853 
0.590268 
0.580154 
0.567899 
0.552456 



Table 6: Comparison of our Monte Carlo results for the magnetization with 
a Fade ' resummation of the series of ref. P^/ . 



3] 



C(t) = ao + A±\t\ 



hg\t\y 



1- - 



3 log(- log |t| 
2 -log|t| 



+ 



ai 



-log|t| 



(40) 

We performed these two (or three) parameter linear fits, first taking 
into account all the data and then systematically eliminating those far- 
thest from the critical point until an acceptable reduced x^ (namely a 
X^ lower than 1) was reached. In all cases, except for the magnetiza- 
tion, with this second step we reached an acceptable C.L. and stopped. 
Notice that in most of the cases such acceptable C.L. could be reached 
keeping all the data, (see the second column of table 8). We also real- 
ized at this stage that the critical amplitudes A± could be obtained in 
a much more efficient way by looking at the internal energy data. We 
shall discuss this point in detail in the next subsection. 

In the next step we added a next-to-leading non-universal correction. 
Among the various possible terms we chose the one giving, in the range 
of values of P of our simulations, the largest contribution. We had to 
resort to this last step only for the magnetization. In this case there are 
two competing corrections. The first one is the 1/log^ term: certainly 
expected due to the presence of the marginal field. But there is also 
second possibility: thanks to the CFT solution of the model we know 
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p 


our MC 


series 


1.130500 


1.2180(9) 


1.2187 


1.120231 


2.1072(23) 


2.1241 


1.117684 


2.5080(26) 


2.5395 


1.115135 


3.0572(36) 


3.1179 


1.112597 


3.8447(51) 


3.9661 


1.110065 


5.0315(72) 


5.3242 


1.107540 


7.053(14) 


7.795 


1.105020 


10.987(19) 


13.647 



Table 7: Comparison of our Monte Carlo results for xi with a Pade' resum- 
mation of the series of ref. ^3^ . 



that in the spectrum there is a subleading magnetic operator to which 
corresponds a new critical index /?' = 3/4. 

The two terms have comparable magnitude, but it turns out that the 
subleading magnetic correction gives a slightly larger contribution in 
the range of interest. So, according to our strategy, we kept only this 
contribution and neglected the 1/log^ one. The resulting scaling func- 
tion is 



m{t) =B\t\-2{- log \t\y 



^ _ 3 log(-log|t|) a]_ 



16 — log|t| — log|t| 



(41) 

Adding also the subleading magnetic term in the fit we found an im- 
pressive lowering of the x^ ■ 

The results of these fits are summarized in table 8. It is important to stress 
that all the quoted errors are statistical. Besides them we also expect system- 
atic errors due to the truncation of the scaling functions to 0{log\t\) or, in the 
case of the magnetization, to the choice of the next-to-leading non-universal 
correction. In sect. 7.2 we shall discuss this problem in more detail. 
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Obs. 


N 


xi 


C.L. 


Amplitude 


ao 


ai 


02 


x+ 


8 


0.93 


46% 


r+ = 0.0223(14) 


0.05(14) 


6.5(4) 




X- 


8 


0.44 


81% 


r_ = 0.00711(10) 


0.02(1) 


-1.24(10) 




6,+ 


8 


0.48 


82% 


/2,+ = 0.192(4) 




1.35(11) 




6,- 


7 


0.44 


82% 


/2,_ = 0.088(4) 




1.06(27) 




m 


8 


0.80 


54% 


B = 1.1621(11) 




-0.220(6) 


-0.144(9) 



Table 8: Results of the fits for susceptibility, magnetization and correlation 
length. In the second column we report the number of data taken into account 
in the fit, in the third column the reduced chi square and in the fourth column 
the confidence level of the fit. The last four columns contain the best fit 
estimates of the parameters of the fit. For the magnetization we also have 
the contribution of a next to leading magnetization operator (see text) 02- 



7.1 The critical amplitudes A± 

The most efficient way to obtain the critical amplitudes A± is to ffi the 
internal energy (for which we have very precise data) instead of the specific 
heat. The bulk value of the energy is known (from duality) to be E{j3c) = | 
however, due to the finite size of the lattice that we simulated we must 
account for possible small deviations from this asymptotic result. We end up 
with the following ffiting function 



E{t) = - + a_i + ao|t| + 3A±|t|3(-lo^ 



>-i 



log|t|) 



^ 3 log 

2 -log|t 



ai 



-log|t| 



2) 

A severe constraint on the results of this ffi is represented by duality which 
implies A^ = A^ . The result of the ffis in the low and high T phases 
are reported in table 9, where it can be seen that the values of Oq, and A± 
extracted in the two phases are indeed compatible within the errors and that 
a_i is, as expected, very small. Combining the two estimates of A^ and A_ 
we extract as our final estimate: A± = 1.30(6). 



Obs. 


N 


xi 


C.L. 


Amplitude 


a_i 


ao 


ai 




7 

7 


0.67 

0.87 


57% 
45% 


A+ = 1.29(5) 
A_ = 1.316(9) 


-0.0001(29) 
-0.0056(5) 


3.0(1.3) 
1.30(18) 


-1.47(48) 
-0.93(6) 



Table 9: Results of the fits for the internal energy. 



7.2 Non- universal corrections 

A crucial ingredient to test the reliability of the above fits is given by the 
magnitude of the non-universal corrections. In the range of (3 values that we 
studied the log of the reduced temperature t takes values which range from 
—4 up to —5. The non-linear contributions listed in tables 8 and 9 must be 
compared with this reference scale. One easily realizes that for all quantities 
these corrections are rather large (they are more or less of the same order of 
magnitude of the universal corrections evaluated in [^] ) and in the particular 
case of the high temperature susceptibility they are very large. This suggests 
that, even if the fits have a very good confidence level, caution is needed in 
assuming the best fit estimates for the amplitudes which could be affected 
by systematic deviations. Notice that there is no hope to control such large 
non-universal contributions by tuning (3 towards the critical temperature. In 
fact it would be necessary to gain at least a factor 10 in log(t) which, as it 
can be easily seen, would imply a huge enhancement of ^. 

We tried to estimate the systematic errors which affect our estimates of 
the critical amplitude with the following method. We repeated the analysis 
discussed above adding in the fitting functions a term of the form l/{log(t)Y 
(which in the range of values of t that we study is the largest among the 
correction to scaling terms that we neglect) with amplitude equal in modulus 
to the Oi amplitudes listed in tables 8 and 9 and with plus and minus sign. 
The two resulting values for the critical amplitude give a (admittedly rough) 
idea of the systematic deviations that we may expect in our estimates. The 
results are collected in table 10. By comparing with the statistical errors 
listed in tables 8 and 9. one can see that the in all cases the systematic 
deviations are larger than the statistical errors. 
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Amplitude 



r+ = 0.0223(40) 

r_ = 0.00711(30) 

/2,+ = 0.192(10) 

/2,- = 0.088(6) 

B = 1.1621(25) 

A± = 1.30(10) 



Table 10: Critical amplitudes with a tentative estimate of the systematic 
errors. 



7.3 Universal Amplitude ratios 

Plugging the values of the critical amplitudes quoted in table 10 in the defi- 



nitions (^, ^ ^, |3^) we find the values for the amplitude ratios reported 
in the second column of table 11. The errors quoted in table 11 have been 
obtained by using the systematic errors quoted in table 10 and discussed in 

sect 7.2. 



8 Comparison with field theory predictions 

It is very interesting to compare our results with the R^, R^,R^, and R^ 
estimates obtained, using the S-matrix form-factors approach to the corre- 
lation functions, in 0. It easy to obtain the remaining four ratios by using 
the following relations: 



R^ 






r; 



R 



R^ 



Ri 



R^ 



(R^y 



Ri 



Rt 

(Rfy 



(43) 



We compare our final estimates and those of in table 11. We immediately 
observe a good overall agreement. This agreement is highly non trivial since, 
as discussed above, in our estimates we had to face large non universal cor- 
rections, while in the predictions of [^ only the lowest states of the spectrum 
were taken into account and some small discrepancies with the exact results 
were expected. If we trust in the overall agreement that we have found we 
immediately see that the major discrepancies between the two sets of data 
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are in the two ratios R^ and -Ri, which could both be consequences of a 
biased estimate of r_. It would be important to understand the reason of 
this discrepancy. In this respect it is worthwhile to notice that xi is the only 
observable for which the non-universal correction has the opposite sign with 
respect to the universal additive one. 



Ratio 


this work 


ref. i 


X 


3.14(70) 


4.013 


R, 


2.19(26) 


1.935 


Ri 


0.68(13) 


0.4539 


R2 


0.44(13) 


0.4845 


Rt 


0.021(5) 


0.0204 


R; 


0.0068(9) 


0.0051 


Rt 


0.220(20) 


0.2052 


R, 


0.100(10) 


0.1060 



Table 11: Comparison between our estimates for the universal amplitude 
ratios and those of ref . ^. 



9 Concluding remarks 

The aim of this paper was to test the recent predictions of for various am- 
plitude ratios in the 4 state Potts model, with the results of a high precision 
Monte Carlo simulation. 

We made four tests on the results of our simulations: 

• Comparison with exact results for small lattices. 



• Comparison with the results of |^, Q at the critical point. 

• Comparison with low temperature series. 

• Agreement with the duality relations. 
All these tests were successfully passed. 
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In looking at the scaling behaviour of our observables, we had to face 
a major problem, due to the presence of a marginal field in the spectrum. 
In performing the analysis we used the recent results of Salas and Sokal on 
the universal additive log log /log correction terms, and found in our fits the 
same behaviour and the same features that they reported in [Q where they 
looked at the finite size corrections at the critical point. 

We found a relatively good overall agreement with the predictions of 
Cardy and Delfino with the exception of two ratios involving the low tem- 
perature susceptibility. It remains an open problem to find a more efficient 
way of dealing, in the analysis of Monte Carlo data, with corrections orig- 
inated by the presence of marginal operators. These contributions are the 
probable main cause of these discrepancies. 
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